Non-quantum Entanglement and a Complete Characterization of pre-Mueller and 

Mueller Matrices in Polarization Optics 



B. Neethi Simon 
Department of Applied Mechanics, IIT Madras, Chennai 600 036 

Sudhavathani Simon 
Department of Computer Science, Women's Christian College, Chennai 600 006 

N. Mukunda 

Centre for High Energy Physics, Indian Institute of Science, Bangalore 560 012 

F. Gori and Massimo Santarsiero 
Dipartimento di Fisica, Universitd Roma Tre and CNISM, Vila della Vasca Navale 84, I-OOI46 Rome 

Riccardo Borghi 

Dipartimento di Electronica Applicata, Universitd Roma Tre and CNISM, Vila della Vasca Navale 84, I-OOI46 Rome 

R. Simon 

The Institute of Mathematical Science, Tharamani, Chennai 600 113 

The Mueller-Stokes formalism which governs conventional polarization optics is formulated for 
plane waves, and thus the only qualification one could demand of a 4 x 4 real matrix M in order 
that it qualifies to be the Mueller matrix of some physical system is that M should map the 
positive cone of Stokes vectors, into itself. In view of growing current interest in the characteri- 
zation of partially coherent partially polarized electromagnetic beams, there is need to extend this 
formalism to such beams wherein the polarization and spatial dependence are generically inseparably 
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physical Mueller matrix. These additional constraints are motivated and fully characterized. 

OCIS codes: 260.5430, 270.5585, 260.2110, 120.5410, 230.5440 



I. INTRODUCTION 



Entanglement is traditionally studied almost exclusively in the context of quantum systems. However, this notion is 
basically kinematic in nature, and so is bound to present itself whenever and wherever the state space of interest is the 
tensor product of two (or more) vector spaces. The vectors of the individual spaces, and hence (tensor) products of 
such vectors, will be expected to possess identifiable physical meaning. Polarization optics of paraxial electromagnetic 
beams happens to have precisely this kind of a setting, and so one should expect entanglement to play a nontrivial 
role in this situation. It turns out that entanglement in this non-quantum setup is not just a matter of academic 
curiosity: we shall show in this paper that this non-quantum entanglement helps in resolving a fundamental issue in 
classical polarization optics. It will appear that this issue could not have been resolved without explicit consideration 
of entanglement. 

A paraxial beam propagating along the positive z-axis is completely determined in terms of the transverse compo- 
nents of the electric field specified throughout a transverse plane z = constant as functions of the transverse variables 
(a;, y) = p. If these components are independent of the transverse coordinates, then the situation corresponds to a 
plane wave propagating along the z-axis. The traditional Mueller-Stokes formalism in terms of Stokes vector S and 
Mueller matrix M, describing respectively the beam and the optical system, presumes essentially this kind of situation 
wherein the spatial degree of freedom can be safely left out of consideration, the focus being on the polarization degree 
of freedom [1,, Q^, ^ . 

Recent years have witnessed an enormous interest in partially polarized partially coherent electromagnetic beams 0, 
S, B 0, B B 0, [m [H, [H, S [3, 0, 0, Hi and hence there is a need to extend the Mueller-Stokes formalism to 



such beams. Given a 4 x 4 real matrix M , it should map $1^?°'), the positive cone of Stokes vectors, into itself in order 
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that it could be the MucUcr matrix of some physical system. Within the conventional formalism,, this seem,s to he 
the only qualification that can be demanded of M. In a partially coherent partially polarized beam, polarization and 
spatial dependence happen to be inseparably intertwined. This inseparability or non-quantum entanglement brings 
in additional constraints that a pre-Mueller matrix M mapping 0(p°') into itself needs to meet in order that it is a 
physically acceptable Mueller matrix. 

The principal purpose of the present work is to motivate these constraints and characterize them fully. The next 
two Sections of the paper act as preparation towards this end. We begin in Section 2 by recounting the Mueller-Stokes 
formalism as it applies to plane waves. This is then extended in Section 3 to paraxial electromagnetic beams, and 
the role of (non-quantum) entanglement between polarization and spatial modulation is rendered transparent. These 
two Sections equip us with all the tools needed to formulate in Section 4 the additional physical constraints on a pre- 
Mueller matrix M, arising as consequence of entanglement. Our final result is formulated in the form of a necessary 
and sufficient condition, and a simple illustrative example is treated in some detail for illustration of the nature of 
these further constraints. And we conclude in Section V with some further remarks. 



II. POLARIZATION OPTICS OF PLANE WAVES 



For a plane wave whose propagation direction is along the (positive) z-axis perpendicular to the {x, y)-plane, the 
X and y components Ei , E2 of the electric field are independent of the transverse- plane coordinates p, and can be 
arranged into a (numerical) column vector 



El 
E2 



(2.1) 



We have suppressed, for convenience, a space-time dependent scalar factor of the form g'^*^^"'^*). While E^E = 
+ I-E2P is (a measure of) the intensity, the ratio 7 = E1/E2 of the (complex) components, which ratio can be 
viewed as a point on the Riemann or Poincare sphere S"^, specifies the polarization state. In particular, the signature 
of the imaginary part of 7 describes the handedness of the (generally elliptic) polarization. 

In presence of fiuctuations E acquires some randomness, and in this case the state of polarization is effectively 
described by the 2x2 coherency or polarization matrix 



$ = {EE"^) 



{EiED {EiED 
{E2EI) {E^E*^) 



(2.2) 



where ( • • • ) denotes ensemble average. The coherency matrix is hermitian, = $, and positive semidefinite, 
V^$y = tr($ W^) > 0, Vy e C^. TMs positlvlty property may be denoted simply $ > 0. Hermiticity and 
positivity are the defining properties of $ : every 2x2 matrix obeying these two conditions is a valid coherency 
matrix, and represents some polarization state. Since $ is a 2 x 2 matrix, the positivity condition takes the simple 
scalar form 



tr$ > 0, 
det $ > . 



(2.3) 



It is clear that the intensity corresponds to tr$. Fully polarized light (pure states) corresponds to det$ = 0, and 
partially polarized or mixed states correspond to dot $ > . 

Typical systems of interest in polarization optics are transversely homogeneous, in the sense that their action is 
independent of the coordinates (a;, y) spanning the transverse plane in which the system lies. If such a system is 
deterministic and acts linearly at the field amplitude level, it is described by a complex 2x2 numerical matrix J 
called the Jones matrix of the system: 



J = 



^11 J12 
J21 J12 



E^ E' 



E[ 
E' 



JE 



^ <^ = {ee'<)^^' = {e'e"') = j<^j'<. 



(2.4) 
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It is clear that Jones systems map pure states (det $ = 0) into pure states. 

Since ^ is hermitian, it can be conveniently described as real linear combination of the four orthogonal hermitian 
matrices tq = l2x2, n = 0-3, ^2 = cri, T3 = 0-2 : 



1 ^ 

$ = 2 XI Sa = tr(ra$); tr r^Tb = 2Sab- 

a=0 



(2.5) 



The reason for choosing the r-matrices, a permuted version of the Pauli matrices rather than the Pauli cr-matrices 
themselves, is to be consistent with the optical convention that the circularly polarised states, the eigenstates of (T2, be 
along the 'third' axis (polar axis) of the Poincare sphere. The intensity equals Sq = tr<i>. The expansion coefficients 
Sa are the components of the Stokes vector S G R*. Note that T3 = — T3 and t* = if a 7^ 3 . 

While hermiticity of $ is equivalent to reality of 5* G R^, the positivity conditions tr$ > 0, det $ > read, 
respectively. 



^0 > 0, 
S'o - S', ^ Si - Si > 0. 



(2.6) 



Thus, permissible polarization states correspond to the positive light cone and its interior (solid cone). Pure states 
live on the surface of this cone. As suggested by this light cone structure, the proper orthochronous Lorentz group 
SO(2>, 1) plays quite an important role in polarization optics 



IE 



Under the action of a deterministic or Jones system J described in (2.4) the elements of the output coherency 
matrix are obviously linear in those of This, in view of the linear relation (2.5) between $ and S*, implies that 
under passage through such a system the output Stokes vector 5" and the input 5* will be linearly related by a 4 x 4 
real matrix M (J) determined by J: 



J: S ^ S' = M{J)S. 



(2.7) 



We may call M{J) the Mueller matrix of the Jones system J. It is known also as a Mueller- Jones matrix to emphasise 
the fact that it is constructed out of a Jones matrix. While $ — (EE^) is a 2 x 2 matrix, the tensor product 
^ = (E ^ E*) is a four-dimcnsional column vector associated with $: 



$ = 



" $0 " 




r (^i^^r) 1 




■$11 - 


$1 




(EiE*) 




«'12 


$2 




{E2E1) 




^'21 


. ^3 . 








. ^'22 . 



(2., 



This idea of going from a pair of indices, each running over 1 and 2, to a single index running over to 3 and vice 
versa can often be used to advantage to associate with any 2x2 matrix K a corresponding column vector K with 
Kq = Kii, Ki = K12, K2 — K21, and K3 = K22- The tensor product J(g) J* is a 4 x 4 matrix: 



J(8) J* 



JuJ* J12J* 
J21J* J22J* 

JllJll JllJ*2 J12J11 J12J12 

J11J21 J11J22 J12J2I J12J22 

J21J11 J21J12 J22J11 J22J12 

.J21J2I J21J22 J22J2I J22J22 



(2.9) 



The transformation (f> 
through 



$' — J<i) jt is thus equivalent to l> $' = J (g) J*l>. Since $ is related to the Stokes vector 



(2.10) 







- 1 








1 ■ 




" $0 " 


Si 




1 








-1 




?i 


S2 
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1 









ISsl 
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—i 







. ^3 . 
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it follows that 

M{J)=A{J®r)A-^, (2.11) 

A being the 4x4 matrix exhibited in Eq. (2.10); this matrix is essentially unitary: A^^ = \A'^ . 

If det J is of unit magnitude, then M{J) computed by this prescription is an element of 5*0(3,1), the proper 
orthochronous group of Lorentz transformations; this was to be expected in view of the two-to-one homomorphism 
between SL{2, C) and 5*0(3, 1). It follows that for any nonsingular J the associated Mueller- Jones matrix M{J) is 
|det J| times an element of 5*0(3, 1). The prescription (2.11), though, applies to singular Jones matrices as well. 



A. Mueller matrices arising from Jones systems 



A nondeterministic system is described directly by a Mueller matrix M : 5^5" = MS and, by definition, such a 
Mueller matrix cannot equal M{J) for any 2x2 (Jones) matrix J. Given a Mueller matrix M, how to test if it is a 
Mueller- Jones matrix for some J or, equivalently, how to test if the system described by M is a deterministic or Jones 
system? It turns out that this question which has received much attention 
solution 



24|. 



19 



21 



22l |23| has a simple and elegant 



We go over in some detail the construction underlying this solution, for it plays a central role in our analysis to 
follow. The sixteen 4x4 hermitian matrices Uab — \Ta ® ^ with a, b independently running over the index set 
{0, 1, 2, 3}, form an orthonormal set or basis in the (vector) space of 4 x 4 matrices; these matrices are unitary and 
self-inverses : 



Uab = Ir^^r; = Ul = U-,\ 
tr (UabUcd) = SacSbd, a, 6, c,(i e {1,2,3,4}. 



(2.12) 



[ Complex conjugation of the second factor of the tensor product Ta ^ is suggested by the construct J ® J* in 
(2.11) ]. Thus every {hermitian) 4x4 matrix can be written uniquely as a {reat) linear combination of {Ta ® t^}. Aa 
an important consequence of this fact we have : 

Proposition 1 : There exists a natural one-to-one correspondence between the set of all 4 x 4 real matrices and the 
set of all 4 X 4 hermitian matrices. 

Indeed, a real matrix M and the associated hermitian matrix Hm are related in this simple manner : 



1 ^ 

Hm = 2 X! ^labTa ® ■ 
a, b=0 

Similarly, a hermitian matrix H and the associated real matrix AIh are related through 

( Mh )ab = ^tr {HTa ® T,*) , fl, 6 = 0, 1, 2, 3 . 
It is clear that these relations are inverses of one another. We write these in more detail for later use: 



(2.13) 



(2.14) 



H 



M 



Moo + Mil 
+ Moi + Mio 

Mo2 + Mi2 

-i{Mo3 + Mi3) 
M20 + M21 

+ i(M3o-f M31) 

M22 + M33 

- i(M23 - M32) 



M02 + M12 

+ i{Mo3 + Mi3) 

Moo - Mil 

- Moi + MiQ 

M22 - M33 

+ i(M23 + Af32) 

M20 - M21 
+ l{M3o-M3i) 



A/20 + Af21 

-iiMso + Msi) 

M22 - M33 

-i(M23 + M32) 

Moo - Mil 
+ Moi ~ Mio 

M02 - M12 

- i{Mo3 - M13) 



M22 + M33 

+ i(A/23 -M32) 

M20 - M21 

-~l{M3o~M3i) 

M02 - Afi2 
+ l{Mo3-Mi3) 

Moo + Afii 
- Moi - Mio 



(2.15) 
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This matrix in identical form was first presented in 2J|. Of the sixteen 4x4 matrices Uab, only U20, U21, U30 and C/31 
have nonzero entries at the 13 location, and this explains the entry M20 — M21 — ii^ho — M31) for (Hm )i3. Written 
in detail, the relation (2.14) has the form 



Mh = t: 



Hqq + Hii 


Hqo — Hii 


Hqi + Hio 


-i{Hoi - 


Hio) 


+ H22 + i?33 


+ H22 — H33 


+ H23 + H32 


- liH23 


— H32 


^^00 + Hii 


Hqo — Hii 


Hqi + Hio 


-tiHoi - 


Hio) 


— H22 — 


— H22 + ^33 


— H23 — H32 


+i{H23 - 


H32) 


H02 + H20 


H02 + H2Q 


Ho3 + H30 


-liHas - 


H3q) 


+ His + H31 


— Hi3 — H31 


+ H12 + H21 


+ l{Hi2 


— H21 


i{HQ2 — H20) 


i{Ho2 — H20) 


i{Ho3 — H30) 


Ho3 + H30 




-i{Hi3-H3i) 


+ i{Hi2 -H21) 


— H12 - 


-H21 



(2.16) 

The entry —{{Hqi — -ffio) + i{H23 — H32) for {Mh )i3 is explained by the fact that the nonzero entries of U13 are at 
the 01, 10, 23 and 32 locations. 

Fundamental to the structure of the Mueller-Stokes formalism of polarization optics is the following result : 
Proposition 2 : Given a 4 x 4 real matrix A/, it is a Mueller- Jones matrix iff the associated hermitian matrix Hm 
is a one-dimensional projection. That is iff Hm = JJ^ for some (complex) four-dimensional column vector J. If 
Hm — JJ^ , then the 2x2 matrix J associated with J is the Jones matrix of the deterministic system represented by 
M. 



In the original work 2J| where this proposition was formulated and proved, this hermitian matrix Hm was actually 
denoted N. However, in the present work we use instead the symbol H to emphasise hermiticity, the most important 
property of this matrix. 

Consider now a transformation which is a convex sum of Jones systems: 



a>'=^PfeJWci>JW^ >o, j2pk 

fe=i 



1. 



(2.17) 



fe=i 



This transformation may be realized by a set of n deterministic or Jones systems J*-^-*, J*-^-*, • • • , j'"-* arranged in 
parallel, with a fraction pk of the light going through the fc'^ Jones system J^''\ and all the transformed beams 
combined (incoherently) at the output. It can also be viewed as a fluctuating system which assumes the Jones form 
jC^) with probability pk ■ In either case, it is clear that the Mueller matrix M of this nondctcrministic system and its 
associated hermitian matrix are corresponding convex sums: 



(2.18) 



fe=i 



fc=i 



It is useful to denote by {p^, J^*^^ } the convex sum or ensemble realization represented by Eq. (2.17), or equivalently 
by Eq. (2.18). Obviously, such an ensemble or convex sum realization always leads to a positive semidefinite Hm, 
and positive semidefinite H alone can be realized as a convex sum of projections. Thus as an immediate, and 
mathematically trivial, consequence of Proposition 2 we have 

Corollary : An optical system described by M is realizable as a convex sum or ensemble {pk, J'"^'^ } of Jones systems 
iff the associated Hm > 0. If Hm > the number of Jones systems, n, needed for such a realization satisfies n > r, 
where r is the rank of Hm- There is no upper limit on n if r > 2. 



This corollary is physically important, and has attracted considerable attention 



25 



26 



27 



28 



31 



32|. 
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B. Pre-Mueller matrices and their classification 



Given a 4 x 4 real matrix Af , Proposition 2 gives the necessary and sufficient condition for M to arise as the Mueller 
matrix of some Jones system J. That still leaves open this more general question: how to ascertain if a given matrix 
M is a Mueller matrix? This question has an interesting history which is surprisingly recent. 

In traditional polarization optics, which is formulated for plane waves and not for beams, the state space IS 
the collection of all Stokes vectors: 

n(p°i) = {S eR"^ \Sq>Q, S'^GS > } , 

G = diag(l, -1, -1, -1) 
S^GS = S^-Sf-Sj-S^. (2.19) 

Now G is the 'Lorentz metric', and thus the state space fi^P"'' is the positive (solid) light cone in ii^. Since a physical 
Mueller matrix should necessarily map states into states, our question reduces to one of effectively characterising real 
linear transformations in ii^ which map the positive (solid) light cone into itself. While SO(3, 1) U G5'0(3, 1), where 
50(3, 1) is the proper orthochronous Lorentz group 50(3, 1), is the answer in the case of onto maps, the more general 
into case was raised in Ref. |333 as a serious issue in polarization optics. This issue was formulated as two simple 
conditions that M has to meet [Eqs. (2.29), (2.31) of Ref. 33] ], corresponding to the demand that the intensity and 
degree of polarization of the output be physical for every input pure state. Further, the measured Mueller matrices 
of Howell |3J| were tested for these conditions and violation was found in excess of 10% , a magnitude considerably 



larger than the tolerance suggested by the reported measurements. Ref. [33| thus concluded that the Howell system 
fails to map the positive light cone fl^P°^'> into itself and this is possibly the first time that a verdict of this kind was 
made on some published Mueller matrices. 

Subsequent progress in respect of this issue was quite rapid. In a significant step forward Givens and Kostinski [s^ 
derived, based on an impressive analysis of the spectrum of GM^GM, what appeared to be a necessary and sufficient 
condition for M to map fi^P"') into itself. They analysed the Howell system based on their own condition and concluded 
that their results were "m coincidence with the negative verdict on the Howell matrix delivered in ^33]". Soon after, 
van der Mee 29| derived a more complete set of necessary and sufhcient conditions for M to map into itself; 

the analysis of van der Mee too was based on the spectrum of GM^ GM . 

Decomposition of a Mueller matrix M in various product forms, to gain insight into t he phys ical effects M 



could have on the input polarization state, has been an activity of considerable interest 36|, l37l . l38l . |39|. The im 



portance of obtaining the canonical or normal forms of Mueller matrices under the double-coset transformation 



M LiMLr, Lg, Lr G 50(3, 1) was motivated in Ref. 20|, and it was shown that the theorem of Givens and Kostin- 



ski implied that the canonical form of every nonsingular (real) matrix M which maps JI^p"') into itself is diagonal; i.e., 
M = LiM^^'^Lr where M^^) = diag(do, di, ^2, da), rfo > > ^2 > Mai and Li,Lr £ 50(3,1). It turned out that 
while the result of van der Mee is essentially complete [i^l, that of Givens and Kostinski is incomplete. This means 
the diagonal form M*^^-' of Ref. ^2^ noted above is not the only canonical form for a nonsingular M mapping fi'^P"'' 
into itself; there exists a non-diagonal canonical form Af , and this is the case that was missed by the 'theorem' of 
Givens and Kostinski quoted above. In a remarkably impressive and detailed study Rao et al. 40|, |41| have further 
explored and completed the analysis of van der Mee, leading to a complete solution to the question of canonical form 
for Mueller matrices, under double-coseting by 50(3, 1) elements, raised in Ref. |20|. 

Since these canonical forms play a key role in our analysis below, we list them here in a concise form. Matrices M 
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which map the state space fi(P°') into itself divide into two major and two minor families : 

Type I : M = LeM^^^Lr, Le, € 50(3, 1) , 
M(^) = diag(do, di, d2, rfs), 



Type II 



do>di>d2> Idsl ; 
M = LeM^'^^Lr, Le, Lr e S0{3,1) , 





" do 


do - di 





■ 








di 





d2 








. 








d3. 



Polarizer : 



do > c^i > 0, y/dodi > ^2 > Msl ; 
M = LtM^P°^'>Lr, Li, Lr G SO{3,l). 
do do 0' 



, do>0; 



Pin Map : 



do do 



0. 

M = LeM^'P'^'^Lr, Le, Lr e SO{3, 1) , 

do 

do 



0. 



- do > . 



(2.20) 



Since elements of SO{3, 1) have unit determinant, it follows that da in the Type-I and Type-II cases is positive, 
negative, or zero according as det M is positive, negative, or zero. The M matrices in the third and fourth families 
are manifestly singular. The third family is a Jones system, the associated H matrix being a pro.iection; indeed, 
M(p°'' corresponds to a Jones matrix J whose only nonvanishing element is Jn = y/2do. Finally, the PinMap family 
is named so because M^p™^ produces a fixed output polarization state independent of the input, the intensity of the 
output being independent of the state of polarization of the input. This may be contrasted with M^^°^^: while the 
output in the case of M(p°') has an input-independent state of polarization, the intensity does depend on the state of 
polarization of the input. 

The matrix M^p'") is not a Jones system, but it is a convex sum of such systems. To see this, note that a perfect 
depolarizer represented by the Mueller matrix 



^(depol) 



10 









(2.21) 



is a convex sum of Jones systems; it can be realized, for instance, as equal mixture of systems with Jones matrices 
Ta, a = 0, 1, 2, 3. That M(p'") is a convex sum of Jones systems follows from M(p'") = M(p°')M('1^p°'). Alternatively, 
it is readily seen that Hf^ipin) is an equal sum of two projections, and hence M^p™^ is an equal mixture of the Jones 
systems 



(2.22) 



While M(p'") is reahzed as convex sum of two Jones systems, M'^'^'^p"') cannot be so realized with less than four Jones 
matrices. This follows from the fact that Hj^^idapoi) is of full rank whereas iJ^/(pi„) is of rank two. 

The classification of canonical forms for M matrices as given in (2.20) is complete in the following sense. 
Proposition 3 : Every M matrix which maps the state space f2(P°') into itself falls uniquely in one of the four families 
described in (2.20). 



" 1 


0' 




'0 


1 " 








7 
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That brings us to the main thesis of the present paper. A 4 x 4 real matrix M will have to map the state space fi^P"'' 
into itself in order that it qualifies to be the Mueller matrix of some physical system. This is certainly a necessary 
condition. And, within the conventional Mueller-Stokes formalism, no conceivable further demand can be imposed on 
M. But the action of the transversely homogeneous system represented by the numerical matrix M can be extended 
from plane waves to paraxial beams; naturally, M will then affect only the polarization degree of freedom and act as 
identity on the (transverse) spatial degrees of freedom. If M indeed represented a physical system, even this extended 
action should map physical states into physical states. It turns out that this trivial looking extension is not all that 
trivial: there are M matrices which appear physical at the level of the (restricted) state space f|(P°i), but fail to be 
physical on the extended state space. Our task in the rest of the paper is to identify precisely those M matrices whose 
action is physical even on the extended state space. Since only those M matrices which pass this further hurdle can 
be called physical Mueller matrices, and pending determination of the precise demand this hurdle places on M , the M 
matrices which map fJ^P"') into itself will be called pre-Mueller matrices. We may thus conclude this Section by saying 
that Eq. (2.20) gives a complete classification of pre-Mueller matrices and their orbit structure under double-coseting 
by elements of 50(3, 1) ; the physical/unphysical divide of pre-Mueller matrices remains to be accomplished. This 
divide will be presented in Section 4 after some further preparation in Section 3. 



III. FROM PLANE WAVES TO BEAMS: THE BCP MATRIX 

We will now go beyond plane waves and consider paraxial electromagnetic beams. The simplest (quasi-) 
monochromatic beam field has, in a transverse plane z = constant described by coordinates {x,y) = p, the form 
E{p) = {Eix + E2y) tp{p), where Ei, E2 are complex constants, and the scalar-valued function ip{p) may be assumed 
to be square-integrable over the transverse plane: ip{p) G L'^{B?). It is clear that the polarization part {Eix H- E2y) 
and the spatial dependence or modulation part ')/'(/') of such a beam are well separated, allowing one to focus at- 
tention on one aspect at a time. When one is interested in only the modulation aspect, the pa rt ( EiX -\- E2y) may 



be suppressed, thus leading to 'scalar optics': this is the domain of traditional Fourier optics [43] ■ [Fourier optics 



for electromagnetic beams requires a more delicate formalism [43| . ] On the other hand, if the spatial part '4'{p) is 
suppressed we are led to the traditional polarization optics (of plane waves) described in the previous Section. 

Beams whose polarization and spatial modulation separate in the above manner will be called elementary beams. It 
is clear that elementary beams remain elementary under the action of transversely homogeneous anisotropic systems 
like waveplates and polarizers. That they remain elementary under the action of isotropic or polarization insensitive 
systems like free propagation and lenses is also clear. 

Now suppose we superpose or add two such elementary beam fields {ax -f by) ili{p) and {cx + dy) x{p)- The result 
is not of the elementary form [ex -\- fy) 4>{p)^ for any e, /, 0(p), unless either (a, b) is proportional to (c, d) so that 
one gets committed to a common polarization, or ^/'(p) and x{p) ^-^e proportional so that one gets committed to a 
fixed spatial mode. In other words, the set of elementary fields is not closed under superposition. 

Since one cannot possibly give up superposition principle in optics, one needs to go beyond the set of elementary 
fields and pay attention to the consequences of inseparability or entanglement of polarization and spatial variation 
(modulation). We are thus led to consider (in a transverse plane) beam fields of the more general form E{p) = 
Ei{p)x + E2{p)y. This form is obviously closed under superposition. We may write E{p) as a generalised Jones 
vector 



Ei{p) 
E2{p) 



E(p) = Ei{p)x + E2ip)y ^ E{p) = 

Ei{p), E2ip)eL^iR^). (3.1) 



The intensity at location p corresponds to |i?i(p)p + |i?2(p)p. This field is of the elementary or separable form iff 
El (p) and E2 (p) are linearly dependent (proportional to one another). Otherwise, polarization and spatial modulation 
are inseparably entangled. 
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The point is that the set of possible beam fields in a transverse plane constitutes the tensor product space C'^®L'^{B?'), 
whereas the set of all elementary fields constitutes just the set product x L'^(R'^) of and L'^{R^). [Recall that 
the tensor product of two vector spaces is the closure of their set product under superposition.] Thus the set product 
X L'^[R^) forms a measure zero subset of the tensor product ® L^(ii^). In other words, in a beam field 
represented by a generic element of ® L^{R'^) polarization and spatial modulation should be expected to be 
entangled : Entanglement is not an exception; it is the rule in (g) L^(i?^), the space of pure states appropriate for 
electromagnetic beams. 

If a beam described by generalised Jones vector E{p) with x, y components Ei{p), E2{p) is passed through an 
x-polarizer, it is not only that the output will be s-polarizcd, it is certain to be in the spatial mode Ei{p) as well. 
Similar conclusion holds if the beam is passed through a y-polarizer. Thus a (transversely homogeneous) polarizer, 
whose action is /j-independent, not only chooses a polarization state, hut acts as a spatial mode selector as well. This 
is true even if Ei{p) and E2{p) are not spatially orthogonal modes. In a similar manner, a spatial mode selector 
insensitive to polarization will end up acting also as a polarization discriminater. This is but one ramification of 
inseparability or entanglement between polarization and spatial modulation. 

Now to handle fiuctuating beams, we pass on to the beam-coherence-polarization (BCP) matrix p') = 

{E{p)E{p'Y), defined as the ensemble average of an outer-product of (generafised) Jones vectors 0, Q. As the 
name suggests, the BCP matrix describes both the coherence and polarization properties of the beam under consid- 
eration. It is a generalization of the numerical coherency matrix of plane waves considered in the previous Section, 
Eq. (2.2), now to the case of beam fields. It can equally well be viewed as a generalization of the mutual coherence 
function of scalar statistical optics to include polarization. We are free to view the BCP matrix either as the 2x2 
matrix of two-point functions 



$(p;p') 



{Ei{p)E,{p'r 
{E2ip)E,ip'r 



{E,{p)E2{p') 

{E2{p)E2{p') 



(3.2) 



or as the associated column vector ^{p;p') of two-point functions: ^{p;p') — {E{p) (g) E[p')*). For our present 
purpose, there is no need to make any finer distinction between the space-time and space-frequency descriptions. 
The two-point functions appearing in the BCP matrix may be viewed either as equal-time coherence functions or 
correlation functions at a particular frequency. 

It is clear from the very definition (3.2) of BCP matrix that this matrix kernel, viewed as an operator from 
(g) L'^{R?') ^ C"^ ® L^{R^), is hermitian and positive semidefinite: 



E 

jk 



<^jk{p;p') = $fcj(p';p)*, i, A: = 1, 2; 
j d^pd^p'E{p)^^{p-p')E{p')>Q, 

I d^pd^p'E,{p)*^,k{p;p')Eu{p') > 0, 

yE{p) e C^^L^iR^). 



(3.3) 



The positivity requirement thus demands that the expectation value of ^(p; p') be nonnegative for every Jones vector 
E{p). Hcrmiticity and positivity are the defining properties of the BCP matrix: every 2x2 matrix of two-point 
functions ^jk{p] p') meeting just these two conditions is a valid BCP matrix of some beam of light. 



We can use the BCP matrix to define, in an obvious manner, the generalised Stokes vector S{p; p') |l4 |: 



1 ^ 

$(p;p') = - Y^Sa{p;p')Ta 

a=0 



^ 5,(p;p') = tr($(p;p')T,). 



(3.4) 



That this is an invertible relation shows that <&(p; p') and S'(p; p') carry identical information: action of an optical 
system on one defines a unique equivalent action on the other. The hermiticity and positivity requirement on the 
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BCP matrix can be easily transcribed into corresponding requirements on S{p; p') . Hermiticity reads 

Sa{p;p')=S^{p';pr, a = 0,1,2,3; 

wtiereas positivity reads 



3 
a=0 



'^pd^p'Sa{p;p')Sa{p';p)>0, 



(3.5) 



(3.6) 



for every Stokes vector S{p;p') arising from Jones vectors of the form E{p) e Z/^(-R^). The signatures Qa 
correspond to the 'Lorentz metric': go = 1, gfo = — 1 for a ^ 1. 

IV. FROM PRE-MUELLER MATRICES TO MUELLER MATRICES: THE ROLE OF 

ENTANGLEMENT 

We now have at our disposal all the tools necessary to determine if a given pre- Mueller matrix is a physical Mueller 
matrix or not. Let us consider the transformation of the generalized Stokes vector S{p;p') and the associated BCP 
matrix ^{p; p') under the action of a transversely homogeneous optical system described by pre-Mueller matrix M. 
We begin our analysis with pre-Mueller matrices of Type-I. 



A. Type-I pre-Mueller matrices 

We will first study pre-Mueller matrices presented in the canonical form M^^) = diag((io, ^i, c^2, ds). Extension of 
the conclusions to pre-Mueller matrices not in the canonical form will turn out to be quite straightforward. In view 
of the system's homogeneity, the action of M^^^ is necessarily independent of p, p', and we have 

M(i) =diag(do, di, d2, ds) : 



r5o(p;p')1 




lS'oip;p')l 




' do Soip; p') ' 


Si{p;p') 




s'Mp') 




di Siip;p') 


S2{p;p') 




S'2{p;p') 
lS's{p;p')\ 




d2S2{p;p') 


lSs{p;p')\ 






.dsSs{p;p'). 



(4.1) 

The elements of the output BCP matrix ^'{p;p') associated with the output Stokes vector S'{p;p') resulting from 
the action of M^^^ on S{p;p'), are easily computed using Eq. (3.4): 

<f>'ii(p;p') = [(do + di)<fn(p;p') + W)-rii)$22(p;p')]/2, 

*22(p;p') = [(do + 'ii)$22(p;p') + (rfo-rfi)*ii(p;p')]/2, 

*i2(P;P') = [(t^2+d3)$12(p;p') + (rf2-rf3)$2l(p;p')]/2, 

^21 (P;p') = [{d2 + d3)^2i{p;p') + {d2-d3)^i2{p;p')]/2. 

(4.2) 

Clearly, a necessary condition for the pre-Mueller matrix M^^^ = diag((io, di, d2, ds) to be a physical Mueller matrix 
is that the output $'(p; p') in Eq. (4.2) be a valid BCP matrix, for every valid input BCP matrix $(p; p'). Hermiticity 
of $'(p; p') is manifest in view of that of $(p; p') and reality of the parameters da- Thus what remains to be checked 
is the positivity of $'(p; p'). While testing positivity of a generic matrix kernel could be a formidable task in general, 
it turns out that this test can be carried out fairly easily in the present case. 

Let us take as the input a special pure state BCP matrix ^^"^ (p; p') = E{p)E{p') \ corresponding to the generalised 
Jones vector E{p) which is an equal superposition of an a;-polarized mode and a y-polarized mode, the two modes 
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being spatially orthogonal: 



E{p) 



= E{p)E{p')'^ 

V'i(p) 



V'2(P) 



(4.3) 



[ipi{p) and ijj2{p) could, for instance, be two distinct Hermite-Gaussian modes.] This means that the entries of 
the input BCP matrix $('')(p;p') have the deterministic form $^°-'(p;p') = ■>pj{p)ijjkip')* ■ A consequence of this 
simple (product) form is that the four entries of the BCP matrix $(°)(p; p') form an orthonormal set of (two-point) 
functions : 



/ 



d^pd^p'^f;\p;p')<i><-^^>{p;pr 



ikOjl- 



(4.4) 



This fact will prove to be of much value in our analysis below. 

To test positivity of the output BCP matrix ^'{p;p'), given in Eq. (4.2) and resulting from input $^°)(p;p') = 
E{p)E{p'y , let us define four (generahzed) Jones vectors: 



E^^\p) = 



Mp) 
±Mp) 



, FW(p) 



V'2(P) 

±Mp) 



(4.5) 



[The input Jones vector E(p) happens to coincide with E^^\p) ]. Expectation values of $'(p; p') for the four Jones 
vectors E'^^\p), F'^^\p) are easily computed using Eqs. (4.2), (4.3), (4.4) and (4.5) in Eq. (3.3). These expectation 
values are {do + di) ± (^2 + ^3) for E'^^\p) and {do - di) ± (^2 - rfs) for F'^^Xp). 

Now positivity of ^'{p;p') requires, as a necessary condition, that these four expectation values be nonnegative, 
and this demand places on the parameters da the constraints 

- rfi - ^2 - ^3 < do, 

-di +d2 + d3 < do; 

di+d2- d3 < do, 

di-d2 + d3 < do. (4.6) 

Violation of any one of these four conditions will render the output $'(p; p') unphysical as BCP matrix. Since 
the input BCP matrix $(°)(p;p') is obviously physical, this will in turn render M^^^ unphysical as Mueller matrix: 
Eq. (4.6) is thus a set of necessary conditions for the pre-Mueller matrix M*-^) to be a Mueller matrix. 

Suppose these four inequalities are met. Can we then conclude that the pre-Mueller matrix M^^^ is a physically 
acceptable Mueller matrix? To answer this question in the affirmative we write in detail the associated hermitian 

matrix Hj^a) = ^J2a ^aTa T* = J2a daUaa ■ 



H. 



M 



do+di 

do — di 

^2 - ^3 

d2 + d3 



^2+^3 

d2 - da 

do — di 

do + di 



(4.7) 



Validity of the four inequalities in Eq. (4.6) implies, fortunately, that this matrix is positive semidefinite. This in 
turn implies that the given diagonal system M^^^ is a convex sum of Jones systems, and therefore takes every BCP 
matrix into a BCP matrix, showing that Eq. (4.6) is sufficient condition for M^^) to be a Mueller matrix. We have 
thus proved 

Proposition 4 : The pre-Mueller matrix M^^^ = diag (do, di, d2, ds) is a Mueller matrix iff the associated hermitian 
matrix Hj^{i) > 0. That is, iff M^^^ can be realized as a convex sum of Jones systems or, equivalently, iff the entries 
of M^^^ respect the inequalities in (4.6). 

Having settled the diagonal case, we now go beyond and consider non-diagonal Type-I pre-Mueller matrices. As 
noted in (2.20), these are necessarily of the general form M = LgM^^'^Lr, where Lg, Lr G SO{3, 1) and M^^^ is 
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diagonal. Wc have already noted that Li, Lr are physical Mueller matrices: indeed, they correspond to deterministic 
systems with respective Jones matrices Ji, Jr € SL{2, C). Thus if M^^^ is a Mueller matrix, then it has a convex sum 
realization {pfe, J^^^ }. This implies that M = LiM^^^'L^ has the convex sum realization {p^, JiJ'^'^^Jr }, and hence 
is a valid Mueller matrix. The converse follows by virtue of the invertibility of J^, Jr, and we have 
Proposition 5 : A Type-I pre-Mueller matrix, which is necessarily of the form M = LfM^^^^L^ with Lg, L^. S 50(3, 1) 
and M*^^) diagonal, is a physical Mueller matrix iff M^^^ is. 



B. Type-II pre-Mueller matrices 



Having fully classified Type-I pre-Mueller matrices into Mueller and non- Mueller matrices, we now turn our attention 
to Type-II pre-Mueller matrices. The analysis turns out to be quite parallel to the one in the previous sub-Section. 

Recall from Section 2 that a Type-II pre-Mueller matrix in its canonical form M^"^^ has only one nonvanishing off- 
diagonal element whose value is fixed by the diagonals, namely Mqi = do — di, where do, di, d2, d^ are the diagonals. 
The action of M^^^ on S{p;p') and <^{p;p') can be computed as before. The (generalised) Stokes vector has this 
simple transformation law: 



M(2): S{p;p') ^ S'{p;p') 



s'o{p;p') 
s'Mp') 

S'^{p;p') 

do So{p; p') + {do - di) 5i(p; p') 
diS,ip;p') 
d2S2ip;p') 
dzS^{p-p') 

The elements of the output BCP matrix $'(p; p') associated with S\p; p') and computed from (3.4) are 

*ii(p;p') = do^ii{p]p'), 

*22(P;P') = di^22{p;p') + {do-di)^ii{p;p'), 

^'i2(p;p') = [(d2 + d3)<i>i2(p;p') + (rf2-d3)$2i(p;p')]/2, 
<&2i (P;p') = [(rf2 + d3)*2i(p;p') + (c!2-ci3)$i2(p;p')]/2. 



(4.E 



(4.9) 



As in the case of M'^^\ the canonical form pre-Mueller matrix M^^) does not couple the pair $ii(p;p'), $22(p;p') 
with $i2(p;p'). «'2i(p;p')- 

Again, a necessary condition for the pre-Mueller matrix M^^^ to be a physically acceptable Mueller matrix is that the 
output ^'{p\ p') in Eq. (4.9) be a valid BCP matrix for every valid input BCP matrix $(p; p'). As in the case of M^^^ 
let us take as input the pure state BCP matrix <I>(°)(p; p') = E{p)E{p') \ with E{p) as described in Eq. (4.3). To test 
positivity of the output BCP matrix $'(p; p'), given in Eq. (4.9) and resulting from input $'^"'(p; p') — E{p)E[p')\ 
we use in place of E^^\p) and i^^^-'(p) slightly modified (generalized) Jones vectors E'^^\p) and F'^^\p) : 



E^'\p) = 



cos 9 tpi (p) 
sin 9 tp2 (p) 



F(^)(p) = 



cos 9 ■ip2 (p) 
sin 9 ipi (p) 



(4.10) 



Expectation values of the output BCP matrix $'(p; p') for these two families of Jones vectors can be computed as 
before. These expectation values are do cos^ 9 + di sin^ ^ + (^2+^3) cos 9 sin 9 for E\p) and (do - di) sin^ 9 + {d2 - 
ds) cos 61 sin 6* for -F^(p). 

Now positivity of $'(p; p') requires, as a necessary condition, that these expectation values be nonnegative for all 
< 9 < TT , and this requirement is seen to be equivalent to the pair of conditions dodi > {d2 + 1^3)^/4, ^2 — 1^3 = 0; 
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these arise respectively from the E^^^ and F^^^ families. We may rewrite these as 

d3 = d2, {d2f<dodi. (4.11) 

This is a pair of necessary conditions for the pre-Mueller matrix M(^) to be a Mueller matrix. The condition 
(c^2)^ < dodi is already part of the definition of M^'^\ and thus d^ = d2 is the new requirement arising from 
consideration of the action of M*^^^ on BCP matrices, i.e., from consideration of entanglement. 

Our next task is to show that these conditions are sufficient as well. To this end we proceed as in the case of M^^^ 
and compute the hermitian matrix Hj^j(2) associated with M'^'^^ : 

do ^{d2+d3y 

rr _ i(d2-rf3) 

^(^2-^3) rfo-rfl 

_^{d2 + d3) . 

(4.12) 

The inequahtics in Eq. (4.11) are precisely the conditions under which Hj^^(2) is positive. This in turn imphes that 
M^^) satisfying (4.11) is a convex sum of Jones systems, and hence is a Mueller matrix. We have thus proved 
Proposition 6 : The pre-Mueller matrix M^^^ is a Mueller matrix iff the associated hermitian matrix H^(2) > 0. 
That is, iff M^^) can be realized as a convex sum of Jones systems or, equivalently, iff the entries of M^^^ respect the 
inequalities in (4.11). 

We can now proceed to consider Type-11 pre-Mueller matrices which are not of the canonical form M^^^ . We know 
from Section 2 that any such matrix has the form M = L^M^'^^Lr, where L^, € 50(3, 1). By considerations similar 
to the ones leading to Proposition 6 in the Type-I case, we arrive at 

Proposition 7 : A type-II pre-Mueller matrix, which is necessarily of the form M = L^M^'^^L^ with L^, € 50(3, 1), 
is a Mueller matrix iff M^^^ is. 

Having completed classification of the pre-Mueller matrices in the Type-I and Type-ll families into physical and 
non-physical ones, we are now left with two minor families to handle. As noted following (2.20), M'p"'' is a Jones 
system. Let J be the Jones matrix representing this system (polariser). In view of the two-to-one homomorphism 
between SL{2,C) and 50(3,1) alluded to earlier, Li, Lr G 50(3, 1) define respective Jones matrices Ji, Jr of 
unit determinant; these Jones matrices are unique except for multiplicative factor ±1 and, as is well known, this 
signature ambiguity is of nontrivial origin. Thus M = LiM^P°^^Lr is a Jones system with Jones matrix ±j£jJr, 
and hence is physical. Similar argument will show that the last family, namely the PinMap family, too has no 
non-physical M matrix. For completeness, we state the situation in respect of these two minor families as the following 

Proposition 8 : Pre-Mueller matrices belonging to the Polarizer and Pin Map families are respectively Jones systems 
and convex sums of Jones systems. Their associated H matrices are positive semidefinite, and all pre-Mueller matrices 
in these two families are physical Mueller matrices. 

C. Complete Characterization of Mueller Matrices 

In the last two sub-Sections we carried out a complete classification of pre-Mueller matrices into physical and non- 
physical ones. Double-coseting under the 50(3, 1) group has played such an important role in this process that we 

capture this role as a separate result. 

Proposition 9 : Given two 4x4 real matrices M and M' which are in the same double-coset orbit under 50(3, 1), 
i.e., M' = L^MLr for some Lc, Lr G 50(3, 1), M' is a convex sum of Jones systems iff M is. And Hm' > iff 
Hm > 0. In other words, M' is a Mueller matrix iff M is. 

Proof: Suppose M has the convex sum realization {pk, J^'^-'}, i.e., M = '^^Pk M{J^''^). Then, clearly, M' has the 
convex sum realization {pk, JiJ^'^^ Jr\- Conversely, if M' has the convex sum realization {p^, J''^'^^}, then M has the 
convex sum realization {p^, ( J^)" V ('^^ ( J^)"!}. 
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Now suppose Hm > 0. This means Hm = Pfe •/^'^^ ^^'^^ ^ Pk > 0. This immediately implies Hm' = 

"^kPk {JiJ'^^^ Jr){JiJ^'^^ Jr) , which proves its positivity. Here {JgJ^'^hJr), as usual, denotes the cohimn vector asso- 
ciated with the 2x2 matrix J^J^''^ Jr. The converse follows from the invertibility of J^, Jr, completing proof of the 
Proposition. 

With this proof of the principal conclusion of this paper is complete. Our main theorem may thus be stated as 
follows. 

Theorem : A 4 x 4 real matrix M is a Mueller matrix iff the associated hermitian matrix Hm > 0. Every physically 
acceptable Mueller matrix is a convex sum of Mueller- Jones matrices. 



D. The role of entanglement : An illustrative example 



We present a simple example to illustrate the kind or restrictions on M matrices brought in by consideration of 
entanglement. Let us restrict attention to M matrices of the special simple three-parameter form 



M 



10 

di 

da 

ds 



(4.13) 



We are obviously in the Type-I situation, but we are not considering here the SO{3, 1) orbit under double-coseting. 
It is clear that M will map Stokes vectors into Stokes vectors if and only if M satisfies the following three conditions: 



-1 < dk < 1, 



fc = 1, 2, 3. 



(4.14) 



In the absence of considerations of entanglement these would have been the only conditions M will need to satisfy. 
Thus the allowed region in the Euclidean space spanned by the parameters {di, d2, ds) would have been the cube 
with vortices at (±1, ±1, ±1). 

Now each one of the four conditions in (4.6) with do = I, arising out of consideration of entanglement, forbids the 
(open) half-space on one side of a plane. Thus the allowed region is the intersection of the four allowed half-spaces. 
This region is clearly the tetrahedron with vertices at (1, 1, 1), (1, —1, —1), (—1, 1, —1) and (—1, —1, 1). 

Thus, it is not the entire cubical region (4.14) permitted by conventional wisdom, but only the solid tetrahedron, 
with one third volume of the cube, that stands the closer scrutiny presented by consideration of entanglement. This 
is illustrated in Fig.l. The region outside the tetrahedron is unphysical and does not correspond to Mueller matrices: 
points in the cube but exterior to the tetrahedron correspond to pre-Mueller matrices which are not physical Mueller 
matrices; those outside the cube correspond to M matrices which are not even pre-Mueller matrices. 



V. CONCLUDING REMARKS 

We conclude with some further observations. In the mathematics literature and in the literature of quantum 
information theory, what we have called pre-Mueller matrices go by the name positive maps, and the subset of pre- 
Mueller matrices which are physically acceptable in the sense of our main theorem corresponds to what are called 
completely positive maps. But we have cndcvoured here to arrive at a physical characterization of Mueller matrices 
entirely within the framework of BCP matrices familiar to the classical optics community, without resorting to the 
mathematical theory of these maps. 

Secondly, a pre-Mueller matrix which fails our main theorem will not produce any unphysical effect acting on the 
coherency matrix of plane waves or on the BCP matrix of elementary (or polarization-modulation separable) beams. 
It follows that BCP matrices which are convex sums of elementary beams will not be able to witness the failure of 
a pre-Mueller matrix M whose associated Hm is not positive semidefinite. Only BCP matrices corresponding to 
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entangled generalized Jones vectors can expose the unphysical nature of a pre-Mueller matrix which violates our main 
theorem. In other words, pre-Mueller matrices cannot be further divided into physical and unphysical subsets without 
consideration of entanglement. 

Finally, ever since it was proved that every Jones system corresponds to a Mueller matrix whose associated H matrix 
is a projection ^2?!, it has been clear that ensembles of Jones systems necessarily correspond to positive semidefinite 




H matrices, and conversely. It has thus been occationally suggested by various authors, beginning with [25|, that 
considerations of Mueller matrices might be restricted to only such ensembles. But it has remained only a suggestion, 
and one without any physical basis, and hence could not set aside as unphysical an experimentally measured Mueller 
matrix whose associated H matrix has a negative eigenvalue, particularly when the reported Mueller system was not 
realized by the experimenter specifically as an ensemble of Jones systems. For instance the symmetric Mueller matrix 
of van Zyl was analysed in Ref. |20|, and was found to be a matrix of Type-I, with canonical- form parameters 
do = 0.9735, di ^ 0.9112, da 0.4640, ds = -0.3838. This clearly violates [the second constraint in] Eq. (4.6) by a 
substantial extent. Equivalently, the eigenvalues oi Hm are 1.0906, 0.8393, 0.4526, —0.3825. 

That Hm in this case is not positive, and hence the van Zyl system is not a convex sum of Jones systems was 
always known. However, one did not have hitherto a physical basis on which this M could be judged as unphysical. 
Non-quantum entanglement has now given us such a physical basis. 
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